The field inside a random distribution of parallel dipoles. 
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We determine the probability distribution for the field inside a random uniform distribution of 
electric or magnetic dipoles. For parallel dipoles, simulations and an analytical derivation show 
that although the average contribution from any spherical shell around the probe position vanishes, 
the Levy stable distribution of the field is symmetric around a non-vanishing field amplitude. In 
addition we show how omission of contributions from a small volume around the probe leads to a 
field distribution with a vanishing mean, which, in the limit of vanishing excluded volume, converges 
to the shifted distribution. 
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The z-component of the field F z at the origin due to 
an electric or magnetic dipole located at r is given by the 
expression 



F, 



= C±((z.h) 



3(f • z)(r ■ n)) , 



(1) 



where C = d/Aire^ for an electric dipole d = dn and 
C = pom/Air for a magnetic dipole m = mh. z and f 
are unit vectors along the z-axis and r, respectively. The 
field at a location within a random uniform distribution 
of many dipoles is a superposition of terms like the one 
in Eq.(l). The field component from a dipole parallel to 
the z-axis located at a distance r and at a direction 9 
with respect to the z-axis is C(l — 3 cos 2 9)/r 3 , and one 
sees that the average of this expression over directions in 
space vanishes for all distances r. It is hence surprising, 
that the field distribution in Fig. ^ obtained by numer- 
ical simulation, is symmetrical around a non-vanishing 
value of the field. We shall prove analytically that the 
distribution is a shifted Lorentzian, shown as the solid 
line in the figure, and that this is the mathematical limit 
of distribution functions which all have vanishing mean 
values but larger and larger variances. 

The fields from electric and magnetic dipoles give rise 
to the most important interactions of neutral matter, 
and they play significant roles in atomic, molecular, and 
many-body physics. In the conclusion we shall list top- 
ics in current quantum gas and quantum information re- 
search where the present analysis may have important 
consequences. 

A typical distance between dipoles with a given density 
p is r = (3/47T/9) 1 / 3 , and a corresponding typical field 
strength is F — Cr^ 3 . For notational convenience, we P\,n{9) 
will rewrite Eq. (JTJ in terms of these typical values as 



9 



F 



d(r,n), 



(2) 



where d(r, n) is the geometrical factor of Eq. 0). 

To derive the distribution function Pn(9) for the field 
component within a randomly distributed collection of 




FIG. 1: The probability distribution for g = F z /Fq, the 
scaled z-component of the field inside a random distribution 
of dipoles aligned along the z-axis. The symbols show the 
result of a numerical simulation (based on 10 6 realizations of 
a system with 50.000 individual dipoles, the uncertainty is 
represented by the line thickness) , the solid curve is the exact 
Lorentzian solution for the probability distribution. 



N dipoles we shall first derive the distribution P\,n{9) = 
P{g\r < iV 1 / 3 ^) for the contribution from a single dipole 
within a sphere of radius N l ' 3 ro. The combined field 
due to N dipoles within the same sphere is distributed 
according to the N-th order convolution product 



Pn(9) 



N 



N 



6 (£2 9i~9) II Pi.Ni.9j) dg 3 



(3) 
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The probability distribution Pi jv is calculated as 



J-V3_, 



'•o 



-3 



d(r,h) 



3r 2 dr 
UrJ 



(4) 

where we explicitly integrate over the radial distribution 
of dipoles, and where (•) denotes the expectation value 
with respect to the direction towards the dipole. The 
expression readily incorporates also an average over pos- 
sibly varying directions n of the individual dipoles to be 
only briefly considered below. By a simple substitution, 



we rewrite Eq. J3J as 



Ng- 



D(Ng), 



(5) 



where D(g) is a geometrical factor which depends only 
on the distribution of d: 



D(g) =(|d(f,n)|y o tf(u- 



d(r,n) 
9 



)du 



(6) 



We observe the simple scaling of the probability distri- 
bution for the field of a single dipolc in a large volume 
holding on average N dipoles: P\ t N{g) — NPi t i(Ng). 

For dipoles parallel to the z-axis, d(r,n) attains the 
value 

d (p) (6) = l-3cos 2 6>, (7) 

from which we find by integration over solid angles that 

j > (p)( g ) = l [2-(2 + ,g) A /T~^ if -2 < <? < 1, 
3\/3 [2 otherwise. 

(8) 
The fact that D^ p \g) assumes a constant value of £>£/ = 
2/3V3 w 0.3849 for |p| > 2 follows from © because |d| 
is bounded by 2 and it provides P\ t N with algebraic tails 
proportional to g~ 2 . This is illustrated in Fig. where 
Pj { is compared to the distribution corresponding to a 
step approximation of D^ p > with the same limiting value: 



D^\g) 



for \g\ > 2D C 
otherwise, 



(9) 



where the position of the edge is determined by the nor- 
malization of Pi t i(g). 




FIG. 2: The probability distribution P^{g) = g' 2 D {v) (g) 
for the field contribution from a single dipole parallel to the 
z-axis. The dotted curve is based on the step approximation 
D' (g), with the same asymptotic values as D^ p '(g). 

Due to the g~ 2 algebraic tails, the distribution Pii 
has a divergent variance and an ill-defined mean value. 
This type of problems is addressed by generalized (Levy) 
statistics, see e.g. |l|, and the form of the bulk distribu- 
tion Poo = limjv^oo Pn can be calculated by the general- 
ized central limit theorem, see e.g. chapter 17 of |2|. We 



will, however, calculate P^ directly as the limit of P/v 
to establish a formalism where the effect of an excluded 
volume can also be obtained. 

The simple scaling relation between Pi^(g) and 
Pi,i(sO allows us to express Pn in terms of P\\ by rewrit- 
ing the convolution Eq. J3J) in Fourier space as 



Pv(<?) = / e 



Jkg 



Pm(#)) TO 



N dk 
2^ : 



(10) 



where r denotes the Fourier transform: f(k) — 
J e~ tkg f(g)dg. Eq. !JTT)|) implies that to determine Pat 
in the limit of N — ► oo we must know the dependence of 
Pi ) i(fc) for small k. We first rewrite Pi i i(fc) as 



P M (fc) = / e-^g- 2 D^(g)dg 



D(g)-DW(gj)dg, 



,-ikg g -2 
where the first term is conveniently rewritten as 

P 2|fe|2>„ 1 _ cog ( t ) \ 



(11) 



l-2|fe|£> 0o I - 



t 2 



-dt 



(12) 



In a small k expansion of the second integral of i|ll|) , the 
Qth or( j er term vanishes since g~ 2 D{g) and g~ 2 D l - e '(g) 
are both normalized. D and D 1 - 6 ' are equal for all \g\ > 2 
and since D^ is even, the 1 st order term yields —ikg c 
with g c defined as 



9c 



.'/" 



gPi A (g)dg, 



(13) 



'/o 



for any go > 2. Collecting the two parts we find that 
Pi,i(fc) = 1 — TT-Doolfcl — ikg c + 0{k 2 ). Insertion of the 
expression for D^ p '(g) leads to the value 



&>> = -[ 3 + V31o| 



y/3-1 

V3+1 



0.1598. 



(14) 



To calculate the limit of P/v for N — » oo, we rewrite 
P)l as log(Pjv(fe)) = NlogP hl {k/N). Since Pi,i(0) = 1 
and log(l-fu) = u + 0(u 2 ), the leading terms of the series 
expansion of (Pi i i(/c) — 1) will dominate in the limit of 
N — > oo, so that 

logCP«,(fc)) = -7rP/oo|fc| - ig c k, (15) 

from which the limiting distribution follows directly: 



PUg) 



i 



irT 2 + (g-g c ) 2 - 



(16) 



This Lorcntzian with a half width of T = itD^ and 
a displacement of g c is in excellent agreement with 
our numerical simulations shown in Fig. ^ The half 
width ttDSo' Fq m 5.065 Cp and central value g c Fq ~ 



0.6692 Cp of the field distribution P{F Z ) are both pro- 
portional to the dipole density p, and their ratio is inde- 
pendent of p. 

The shift of the most probable value with respect to 
zero is surprising when one considers the vanishing mean 
contribution from any spherical shell around the origin, 
but it is less surprising when one observes the probability 
distribution for the single dipole contribution, shown in 
Figure 2. This distribution is indeed suggestive of a shift, 
but its mean is ill-defined, and (|13|1 provides the proper 
procedure to obtain g c from Pi.i(g). 

For completeness we note that, in the case of randomly 
oriented dipoles, the factor d(r, n) is given by 

$- r ' = sin 9 1 sin 9 2 sin <f> — 2 cos d\ cos 62, (17) 

where 0\ is the direction of f, 62 is the angle between 
r and n, and </> represents the rotations of n around 
r. Integration over these angles with the appropriate 
probability measure sin(0i) sin(0 2 )d0id02<#/87r yields an 
even function D^ r '(g) with the asymptotic limit Doo = 



sinh (vo) « 0.3450, implying a Lorentzian dis- 



4 T 24 

(r) 

tribution with a half width T = irDoo ~ 1.083 centered 
at zero field. This is in agreement with work by Stoneham 
[3j, who considered a variety of line broadening mecha- 
nisms in solids, and identified a Lorentzian line as the 
result of interaction of a single molecule with disloca- 
tion dipoles. More recently |4j Lorentzian line shapes 
were measured for molecules embedded in low tempera- 
ture glass with a low-density distribution of dynamical 
defects. These results were interpreted in terms of Levy 
stable distributions. 

If the algebraic tails of Pi : i are truncated by some 
mechanism, the distributions have finite variance, and 
our naive estimates of mean values will be valid due to 
the central limit theorem. To investigate whether such a 
truncation entirely removes the more spectacular effect 
identified above, we shall compute the field distribution 
in the case where we will not allow any dipoles inside 
an excluded volume in the form of a sphere of volume 
e/p centered at the origin. Note that e is the average 
number of dipoles that would have been found in the ex- 
cluded volume. The symbols in Fig. [3] show the results 
of simulations performed with dipoles put uniformly at 
random around the origin but outside such excluded vol- 
umes, and as we reduce the excluded volume we observe 
that the probability distributions converge towards the 
shifted Lorentzian. The generalized central limit theorem 
which applies for e = and e — > 00 deals with the conver- 
gence of the distribution function for a sum of more and 
more random variables which all have the same individ- 
ual distribution after a suitable rescaling. Such rescaling 
is not possible for intermediate values of e, which thus 
require a direct calculation of Poo(e,g). 

We consider the field contribution from a single dipole 
placed at random in a spherical shell with outer radius 
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FIG. 3: Distribution of g for the case of dipoles parallel to 
the z-axis when an excluded volume of size e/p is introduced. 
Solid lines show the distribution Poo (e, g) given by Ijl8^ , for 
e = 0, 0.4, 1, and 2 in order of increasing maximum densities. 
Data markers are the result of numerical simulations. The 
distributions have vanishing mean for all values of e > 0, yet 
they approach the shifted Lorentzian, corresponding to e = 0. 



(N + e) 1 ' 3 ^ and inner radius e 1 ' 3 ^. Parametrizing the 
radius by x = (r/ro) 3 , the mean number of atoms pop- 
ulating the sphere with radius r, we have by Bayes rule 
and the additivity of probabilities of disjoint events that 
P(e < x < N + e)P{g\e < x < N + e) = P(x < 
N + e)P(g\x < N + e) - P(x < e)P(g\x < e), where 
the first factors are simply the probabilities that a single 
particle is found in the specified regions of space, and, 
e.g., P(x < {N + e))/P{x < e) = (N + e)/e. Taking 
the Fourier transform with respect to g and noting that 
P(k\x < Xi) — Pi,i(k/xi), we obtain the following rela- 
tion between the Fourier transformed probabilities 

N P(*|e < x < N + e) = {N + e) PiAwrJ ~ ^1,1(7)- 

We are interested in the probability PN{e,g) that the 
contributions from N dipoles, all having e < x < N + e, 
add up to the value g. Performing the convolution in 
Fourier space we find that logFV(e, k) — ^VlogP(fc|e < 

x < N + e), and for N — > 00 we have 



logPo^e, k) = logP^fc) - e(Pi,i(|) - 1). 



(18) 



As shown by Fig. [3] this expression, which can be evalu- 
ated numerically, is in excellent agreement with numeri- 
cal simulations. 

I Pi 1 I < 1, and, by (|15|1 . the term log Poo (fc) will dom- 
inate Eq. I|18fl for k > e, in agreement with our expec- 
tation that P 00 (e,<7) should approach P 00 (5 r ) for e — > 0. 
To consider the limit of e — * 00 we continue the series 
expansion of P^i to find logP^ (e, fc) = — 2/5e _1 fc 2 — 



4i/105. 



-2;„3 



eO((fc/e) 4 ). Since the leading term of 



this expansion will dominate for e ^> 1, we conclude that 
Poo(e, g) asymptotically approaches a Gaussian distribu- 
tion with variance Var(<7) = 4/5 e _1 . 



In summary, we have identified a shifted Lorentzian 
distribution as the probability distribution for the to- 
tal field inside a random distribution of dipoles, and 
we have identified a family of distributions for the case 
where dipoles are not permitted inside an excluded vol- 
ume around the origin. These distributions have vanish- 
ing mean, and they converge to Gaussian distributions 
in the limit of large excluded volumes and towards the 
Lorentzian in the case of small excluded volumes. It 
is not an inconsistency of our results that the shifted 
Lorentzian is approached by distributions with vanishing 
mean: a Lorentzian can be ascribed any mean value de- 
pending on how the upper and lower limits are taken in 
the integral over the distribution. There is in fact rea- 
son to emphasize that the common procedure of fitting a 
spectrum to a Lorentzian may be quite misleading if one 
tries to interpret a frequency shift as the mean value of 
a possible perturbation of the energy of the system. 

The emergence of a non- vanishing field as the most 
likely result of the surroundings of any given dipole, could 
have consequences for material properties. The proper- 
ties of a conventional ferromagnet are controlled by an in- 
terplay of Coulomb forces and the Pauli exclusion princi- 
ple for electrons, which may be conveniently represented 
by a spin-spin interaction term, but in novel materials, 
such as recently produced carbon-nanofoams [3] , the ac- 
tual interaction between separated magnetic dipoles may 
be an important ingredient in the understanding of their 
collective properties. 

Heteronuclear molecules with permanent electric 
dipole moments have been trapped |a| and experiments 
are planned with atomic species with particularly high 
magnetic dipole moments pi Bill to study polar degener- 
ate gas es and new kinds of order and collective dynamics 
[1(1 111! ll2J Il3 | . Mean-field approaches in dipolar degen- 
erate quantum gases may of course be questionable if the 
mean field itself is not well defined. Our work suggest 
that a critical examination of this issue is necessary. 

Highly excited Rydberg atoms in electric and magnetic 
fields interact strongly |14j , and fast quantum computing 
[LJ and single photon generating devices [lfi| have been 
suggested based on the energy shifts in atoms caused by 
the excitation of nearby atoms. Rare-earth ions in crys- 
tals have excited states with permanent electric dipole 
moments, and proposals exist for quantum computing 
within such a system which are also based on large |l7| 
or small (18| shifts in absorption frequency of target ions 
caused by excitation of a nearby control-ion. In the rare- 
earth system Lorentzian broadening of spectrally hole 
burnt structures has been observed when ions at different 
frequencies are excited J19j , and we imagine that this can 
be an ideal system to study the broadening and the shift 
systematically, as the density of perturbing dipoles can 
be varied by the exciting laser system. 



Extension of the analysis, e.g., to time-dependent fields 
and to higher order multipole fields seems very inter- 
esting. Von Neuman and Chandrasekhar considered the 
fluctuating gravitational forces in a stellar medium, see 
|20|. As a curiosity we note that the time derivative 
of these forces at any given time behaves like a sum of 
dipole fields, and hence a massive object moving through 
a static random mass distribution may experience a force 
with a time derivative given by the shifted Lorentzian. 
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was funded by the European Union IST-FET programme 
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